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Abstract. We express quantum computations (with measurements) us- 
ing the arrow calculus extended with monadic constructions. This frame- 
work expresses quantum programming using well-understood and famil- 
iar classical patterns for programming in the presence of computational 
effects. In addition, the five laws of the arrow calculus provide a conve- 
nient framework for equational reasoning about quantum computations 
that include measurements. 



1 Introduction 

Quantum computation [1] can be understood as a transformation of information 
encoded in the state of a quantum physical system. Its basic idea is to encode data 
using quantum bits (qubits). Differently from the classical bit, a qubit can be 
in a superposition of basic states leading to "quantum parallelism." This form of 
parallelism is due to the non-local wave character of quantum information and is 
qualitatively different from the classical notion of parallelism. This characteristic 
of quantum computation can greatly increase the processing speed of algorithms. 
However, quantum data types are computationally very powerful not only due 
to superposition. There are other odd properties like measurement, in which the 
observed part of the quantum state and every other part that is entangled with 
it immediately lose their wave character. 

These interesting properties have led to the development of very efffcient 
quantum algorithms, Hke Shor's quantum algorithm for factorizing integers [2], 
and Grover's quantum search on databases ^l- Another important theme is the 
development of quantum cryptographic techniques p] . 

Since these discoveries, much research has been done on quantum computa- 
tion. Summarizing the field of research we can classify it according three main 
areas: i) physical implementations of quantum computers, ii) development of 
new quantum algorithms; and iii) design of quantum programming languages. 

This work is about the design of a quantum programming language, and 
consequently about a high-level, structured and well-defined way to develop new 
quantum algorithms and to reason about them. 



We have been working on semantic models for quantum programming. In 
previous work [5] we established that general quantum computations (including 
measurements) are an instance of the category-theoretic concept of arrows [6], a 
generalization of monads \J\ and idioms [8\ . Translating this insight to a practical 
programming paradigm has been difficult however. On one hand, directly using 
arrows is highly non-intuitive, requiring programming in the so-called "point- 
free" style where intermediate computations are manipulated without giving 
them names. Furthermore reasoning about arrow programs uses nine, somewhat 
idiosyncratic laws. 

In recent work, Lindley et. al. [9\ present the arrow calculus, which is a 
more friendly version of the original presentation of arrows. The arrow calculus 
augment the simply typed lambda calculus with four constructs satisfying five 
laws. Two of these constructs resemble function abstraction and application, and 
satisfy familiar beta and eta laws. The remaining two constructs resemble the 
unit and bind of a monad, and satisfy left unit, right unit, and associativity laws. 
Basically, using the arrow calculus we can understand arrows through classic 
well-known patterns. 

In this work we propose to express quantum computations using the arrow 
calculus axtended with monadic constructions. We show that quantum program- 
ming can be expressed using well-understood and familiar classical patterns for 
programming in the presence of computational effects. Interestingly, the five laws 
of the arrow calculus provide a convenient framework for equational reasoning 
about quantum computations (including measurements). 

This work is organized as follows. The next two sections review the back- 
ground material on modeling quantum computation using classical arrows. Sec- 
tion |4] presents the arrow calculus. We show the quantum arrow calculus in Sec- 
tion m We express some traditional examples of quantum computations using 
the quantum calculus. Additionally, we illustrate how we can use the calcu- 
lus to reason about quantum programs. Section [6] concludes with a discussion of 
some related works. Finally, Appendix [A] presents the constructs of simply- typed 
lambda calculus. Appendix [B] gives an extension of the simply-typed lambda 
calculus with monadic constructions, and Appendix [C] reviews general quantum 
computations. 

2 Classic Arrows 

The simply-typed lambda calculus is an appropriate model of pure functional 
programming (see Appendix [Al . The standard way to model programming in 
the presence of effects is to use monads pO] (see Appendix [B]) . Arrows, like 
monads, are used to elegantly program notions of computations in a pure func- 
tional setting. But unlike the situation with monads, which wrap the results of 
computations, arrows wrap the computations themselves. 

From a programming point of view, classic arrows extend the simply-typed 
lambda calculus with one type and three constants satisfying nine laws (see 
Figure 1). The type B denotes a computation that accepts a value of type 
A and returns a value of type B, possibly performing some side effects. The 



three constants are: arr, which promotes a function to a pure arrow with no side 
effects; 3^, which composes two arrows; and first, which extends an arrow to 
act on the first component of a pair leaving the second component unchanged. 

To understand the nine equations, we use some auxihary functions. The func- 
tion second, is Uke first, but acts on the second component of a pair, and fMig, 
appHes arrow / and g to the same argument and then pairs the results. 



Types 

arr :: {A^ B) ^ {A--^ B) 

:; [A-^ B) ^ [B C) ^ [A-^ C) 
first :: {A-^ B) ^ (A X C B x C) 

Definitions 
second : {A B) ~* (C x A C x B) 
second — Xf.arr swap ^ first f arr swap 
{m) : {C A) ^ [C B) ^ {C Ax B) 
(ffis) — Xf.Xg.arr sup ^ first f ^ second g 

Equations 
(-^i) arr id f = f 

(-^2) / ^ arr id = f 



3 Quantum Arrows 

Quantum computation is generally expressed in the framework of a Hilbert space 
(see Appendix \C\ for a short review of that model) . As expressive and as con- 
venient is this framework for mathematical reasoning, it is not easily amenable 
to familiar programming techniques and abstractions. In recent work [5] how- 
ever, we established that this general model of quantum computations (including 
measurements) can be structured using the category-theoretic concept of arrows. 
Figure 2 explains the main ideas which we elaborate on in the remainder of this 
section. 

In the figure, we have added type definitions (i.e, type synonyms) for conve- 
nience. Type Vec A means that a vector is a function mapping elements from 
a vector space orthonormal basis to complex numbers (i.e., to their probabil- 
ity amplitudes). Type Lin represents a linear operator (e.g, a unitary matrix) 
mapping a vector of type ^ to a vector of type B. Note that if we uncurry the 
arguments A and B, it turns exactly into a square matrix (i.e, Vec {A,B)). 



Fig. 1. Classic Arrows 
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Type Dens A stands for density matrices and it is straight to build from Vec. 
Type Super A B means a superoperator mapping a density matrix of type A 
to a density matrix of type B. This type can be understood by interpreting it in 
the same style as Lin. 



Fig. 2. Quantum Arrows 



Type Definitions 

type Vec A = A C 
type Lin AB ^ Vec B 
type Dens A = Vec {A, A) 
type Super AB ^ {A,A) ^ Dens B 
Syntax 

Types A, B, C ::= ... Vec A \ Lin A \ Dens A \ Super A B 

Terms L,M,N ::= ... | return \ ^ | arr \ \ first 
Monadic Definitions 

return : A Vec A 

return ab — \f a —— b then 1.0 else 0.0 
: Vec A^ {A^ Vec B) Vec B 

va / — Xb. 1 {""1 ")(/ t[ b) 
Auxiliary Definitions 

fun2lin : {A ^ B) ^ Lin A B 

fun2lin / = A a.return (/ a) 

((*)) : Vec A Vec B Vec (A, B) 

Vi{*)v2 ~ A (a,b).vi a * V2 b 
Arrow Types and Definitions 

arr : {A ^ B) ^ Super A B 

arr f = fun2lin (A (6i, 62) ^ (/ ,f b2)) 

:: (Super A B) ^ (Super B C) ^ (SuperA C) 

f g = X b.{f b:^g) 

first :: (Super A B) ^ (Super {A x C) {B x C)) 

first f ((61, di), (62, d2)) = permute ((/(&i, b2)){*) return (di, ^2)) 
where permute v ((&i, 62), (di, 0(2)) = ((&i, di), (&2, d2)) 



We have defined in our previous work [5] the arrow operations for quantum 
computations into two levels. First we have proved that pure quantum states (i.e, 
vector states) are an instance of the concept of monads [7|. The definitions of 
the monadic functions are shown in Figure 2. The function return specifies how 
to construct vectors and defines the behavior of an appHcation of matrix to 
a vector. Moreover we have used the auxiliary functions fun2lin, which converts 
a classical (reversible) function to a linear operator, and (*) which is the usual 
tensor product in vector spaces. 

The function arr constructs a quantum superoperator from a pure function 
by applying the function to both vector and its dual. The composition of arrows 
just composes two superoperators using the monadic bind. The function first 
applies the superoperator / to the first component (and its dual) and leaves the 
second component unchanged. 



We have proved in our previous work that this superoperator instance of 
arrows satisfy the required nine equations [5]. 



4 The Arrow Calculus 



In this section we present the arrow calculus [9] and show the translation of the 
calculus to classic arrows (described in Section [2]) and vice versa. The translation 
is important because it essentially corresponds to the denotational semantic 
function for the quantum version of the arrow calculus. The material of this 
section closely follows the original presentation in |9j. 



4.1 The Calculus 

The arrow calculus as shown in Figure 3 extends the core lambda calculus with 
four constructs satisfying five laws. Type A ^ B denotes a computation that 



Fig. 3. Arrow Calculus 



Syntax 

Types A,B,C 
Terms L, M, N 
Commands P,Q,R 
Types 

r-x: Ah QIB 
r h X'x.Q : A-^ B 

r,Ah M : A 



... \A-^B 
... I X'x.Q 

L • P I [M] I let X = P in Q 



r\- L-.A-^B r;A\- A4 -.A 

r;Ah L» M\B 
r;AhP\A r;A,x : Ah Q\B 
P; Zi h let I = P in Q\B 

= Q[x ~ M] 
= L 

Q = Q[x ~ M] 

= P 



r-Ah [M\\A 
Laws 

(Z?-^) {X'x.Q) 

(r^-) X'x.{L.[x]) 

(left) let a; = [M] in C 

(right) let a; = P in [x] 

(assoc) let y — (let a; = P in Q) in P = let a; = P in (let y — Q tn R) 



accepts a value of type A and returns a value of type B, possibly performing 
some side effects. 

There are two syntactic categories. Terms are ranged over by L, M, N, and 
commands are ranged over by P,Q, R. In addition to the terms of the core 
lambda calculus, there is one new term form: arrow abstraction X'x.Q. There 
are three command forms: arrow application L • M, arrow unit [M] (which 
resembles unit in a monad), and arrow bind let x — P in Q (which resembles 
bind in a monad). 

In addition to the term typing judgment F h M : A there is also a command 
typing judgment r;A\- P\A. An important feature of the arrow calculus is that 



the command type judgment has two environments, F and A, where variables 
in r come from ordinary lambda abstractions Xx.N, while variables in A come 
from arrow abstraction X*x.Q. 

Arrow abstraction converts a command into a term. Arrow abstraction closely 
resembles function abstraction, save that the body Q is a command (rather than 
a term) and the bound variable x goes into the second environment (separated 
from the first by a semicolon). 

Conversely, arrow application, L»M\B embeds a term into a command. Ar- 
row application closely resembles function appHcation. The arrow to be applied 
is denoted by a term, not a command; this is because there is no way to apply 
an arrow that is itself yielded. This is why there are two different environments, 
r and A: variables in F may denote arrows that are applied to arguments, but 
variables in A may not. 

Arrow unit, [M]!A, promotes a term to a command. Note that in the hy- 
pothesis there is a term judgment with one environment (i.e, there is a comma 
between F and A), while in the conclusion there is a command judgment with 
two environments (i.e, there is a semicolon between F and A). 

Lastly, using let, the value returned by a command may be bound. 

Arrow abstraction and application satisfy beta and eta laws, {fS^) and {if^), 
while arrow unit and bind satisfy left unit, right unit, and associativity laws, 
(left) , (right) , and (assoc) . The beta law equates the application of an abstraction 
to a bind; substitution is not part of beta, but instead appears in the left unit 
law. The (assoc) law has the usual side condition, that x is not free in R. 

4.2 Translation 

The translation from the arrow calculus to classic arrows, shown below, gives a 
denotational semantics for the arrow calculus. 

[X'x.Q] = [Q]^ 

[L • M]a = arr{XA.[M]) [L] 

[[M]]a = arr{XA.[M]) 

[let x = P\nQ]A = {arr id Ml [P]a) :» [Q]a,x 

An arrow calculus term judgment F \- M : A maps into a classic arrow judgment 
F h [M] : A, while an arrow calculus command judgment F\ AV- PI A maps into a 

classic arrow judgment F h [P]a : A-^ A. Hence, the denotation of a command 

is an arrow, with arguments corresponding to the environment A and result of 

type A. 

We omitted the translation of the constructs of core lambda calculus as they 
are straightforward homomorphisms. The translation of the arrow abstraction 
X'x.Q just undoes the abstraction and call the interpretation of Q using x. 
Application L • P translates to [M] translates to arr and let a; = P in Q 
translates to pairing Ms (to extend the environment with P) and composition 
»>(to then apply Q). 

The inverse translation, from classic arrows to the arrow calculus is defined 

as: 



[arr]-i = Xf.X'x.[f x] 

K:^)]-! =Xf.Xg.X'x.g.if.x) 

[first]^^ = Xf.X*z.\et a; = / • fst z in [{x, snd z)] 

Again we omitted the translation of the constructs of core lambda calculus as 
they are straightforward homomorphisms. Each of the three constants from clas- 
sic arrows translates to an appropriate term in the arrow calculus. 



5 The Arrow Calculus as a Quantum Programming 
Language 

In this section we discuss how the arrow calculus can be used as a quantum 
programming language. 

We start by showing quantum programs using the standard quantum circuit 
notation. The lines carry quantum bits. The values flow from left to right in steps 
corresponding to the alignment of the boxes which represent quantum gates. 
Gates connected via bullets to another wire are called controlled operations, that 
is, the wire with the bullet conditionally controls the appHcation of the gate. The 
circuit in Figure 4 represents a quantum program for the Toffoli gate. Using the 



Fig. 4. Circuit for the ToffoH gate 
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classic arrows approach for quantum programming presented in Section [3] and 
using the type of booleans, Bool, as the orthonormal basis for the qubit, this 
program would be codded as follows: 

toffoli :: Super (Bool, Bool, Bool) (Bool, Bool, Bool) 

toffoli = arr (A(ao, 6o, co) ^ (co, (ao, 5o))) 

{first H arr (A(ci, (oo,6o)) ((5o, ci), ao))) 
{first cV » arr (A((6i, ca), oq) ^ {{ao , h) , C2))) ^ 
{first cNot arr (A((ai , 62), C2) ^ ((62, C2), ai))) ... 

As already noted by Paterson [TT] this notation is cumbersome for programming. 
This is a "point-free" notation, rather different from the usual way of writing 
functional programs, with A and let. Paterson introduced syntactic sugar for 
arrows, which we have used in our previous work [5]. However, the notation 
simply abbreviates terms built from the three constants, and there is no claim 
about reasoning with arrows. Using the quantum arrow calculus presented in 
Figure m this program would be Hke: 



toffoli :: Super (Bool, Bool, Bool) (Bool, Bool, Bool) 
toffoli = A*.(x, y, z).\et z' = H • z\r\ 
\et{v',z")=cV*{v,z')m 
let (x', y") = cNot • [x, y')in . . . 

This style is more convenient and elegant as it is very similar to the usual 
familiar classical functional programming and is amenable to formal reasoning in 
a convenient way. Consider, for instance, the program which applies the quantum 
not gate twice. That is obviously equivalent to identity. To do such a simple proof 
using the classic arrows we need to learn how to use the nine arrow laws and 
also to recover the definitions of the functions arr, and first for quantum 
computations presented in Figure 2. 

The action of the quantum not gate, QNot, is to swap the amplitude proba- 
bilities of the qubit. For instance, QNot applied to |0) returns |1), and vice versa. 
But QNot applied to a|0) + I3\\) returns a|l) + ^|0). 

Given the classical definition of not as follows: 

not = Ax. if X == True then False else True : Bool Bool 

Using the arrow calculus, the QNot would be written as: 

QNot = A'y.[not y] : Super Bool Bool. 

Then, the program which applies the QNot twice, would be: 

r [- A* a;. let w — (A'z.[not z]) • x\n (A*y.[not y]) • w 

Again the syntax, with arrow abstraction and appHcation, resembles lambda 
calculus. Now we can use the intuitive arrow calculus laws (from Figure 3) to 
prove the obvious equivalence of this program with identity. The proof follows 
the same style of the proofs in classical functional programming. 

A*a;.let w = (A'z.[not z]) • x in (A'y.[not y]) • w 
A'x.let w — [not x] in (A*T/.[not y]) • w 
X'x.{X'y.[not y]) • (not x) 
A*a;.[not(not x)] 
X'x.[x] 

It is interesting to note that we have two ways for defining superoperators. 
The first way is going directly from classical functions to superoperators as we 
did above for not, using the default definition of arr. The other way is going 
from the monadic pure quantum functions to superoperators. As monads are a 
special case of arrows [6] there is always a translation from monadic functions 
to arrows. Hence, any Lin A B is a special case of Super A B. 

Hence, we construct the quantum arrow calculus in Figure 5 in three levels. 
First we inherit all the constructions from simply-typed lambda calculus with 
the type of booleans and with classical let and if (see Appendix [Aj . Then we 



_(left) 

= (/3-") 
def.not 



Fig. 5. Quantum Arrow Calculus 



Syntax 

Types A, B, C 
Terms L,M,N 
Commands P,Q,R 
Monad Types 



= . . . I Bool I Dens A \ Vec A \ Super A B 
= [T] I let X = M in N\ X'x.Q | + | - 
= L» P \ [M]\\etx = P \n Q \ meas | trL 



r\- M -.A Ph M : Vec A r,x : Ah N : Vec B 



r h [M] ■ Vec A r\-\etx = Mm N : Vec B 

Ph M,N : Vec A Ph M,N : Vec A 



P h M+N ■ Vec A r h Af-Af : Vec ^ 
Arrow Types 

P;x: A\-Q\ Dens B P \- L : Super A B T; A h M : A 



rh A'x.Q : Super A B T; A h L • A/! Dens B 

r,AI-A/:A T; Z\ h P! Dens A T; A, x : A h Q! Dens B 



T; A h [A/]! Dens A B; zi h let a; = P in Q! Dens B 



B; a; : A I- meas ! Dens (A, A) B; a; : (A, B) h trL ! Dens B 



add the monadic unit, [], to build pure vectors (over booleans), let to sequence 
computations with vectors, and plus and minus to add and subtract vectors (the 
monadic calculus [7j with its laws is presented in Appendix [BJ . Finally, we add 
the constructions of the arrow calculus. The appeal of using the arrows approach 
is because we can express measurement operations (i.e, extract classical infor- 
mation from the quantum system) inside the formalism. Therefore, we have two 
computations for measurements on mixed states, meas and trL. The computa- 
tion meas returns a classical value and a post-measurement state of the quantum 
system. The computation trL traces out or projects part of the quantum state 
(the denotation of these operations is provided in Appendix [P]) . 

To exemplify the use of the monadic constructions, consider, for example, 
the hadamard quantum gate, which is the source of superpositions. For instance, 
hadamard appHed to |0) returns |0) + |1}, and appHed to |1) returns |0) — |1). 
But, hadamard applied to |0) + |1) returns |0), as it is a reversible gate. To define 
this program in the quantum arrow calculus, we just need to define its work for 
the basic values, |0) and as follows: 

hadamard = Ax. if x =— True then [False] — [True] 

else [False] + [True] : Lin Bool Bool 



Then, the superoperator would be: 



Had = A*y.[hadamard y] : Super Bool Bool 

Another interesting class of operations are the so-called quantum controlled 
operations. For instance, the controlled not, Cnot, receives two qubits and applies 
a not operation on the second qubit depending on the value of the first qubit. 
Again, we just need to define it for the basic quantum values: 

cnot = \{x. y).if X then [(a;, not y)] 

else : Lin (Bool, Bool) (Bool, Bool) 

Again, the superoperator of type Super (Bool, Bool) (Bool, Bool) would be 
Cnot = A*(x, ?;).[cnot {x, y)]. 

The motivation of using superoperators is that we can express measurement 
operations inside of the formalism. One classical example of quantum algorithm 
which requires a measurement operation is the quantum teleportation [4]. It 
allows the transmission of a qubit to a partner with whom is shared an entangled 
pair. Below we define the two partners of a teleportation algorithm. 

Alice : Super (Bool, Bool) (Bool, Bool) 
Alice = A'(a;, y). let (.t', y') = Cnot • [x, y) in 
let q = (Had • x' , y') in 
let [q' , v) = meas •gin trL • (g, v) 

Bob : Super (Bool, Bool, Bool) Bool 

Bob = X'{x, y, z). let (z', x') = Cnot • (z, x) in 

let (y', a;") = (Cz . (y, x')) in trL . ((y', z'), 

6 Conclusion 

We have presented a lambda calculus for general quantum programming that 
builds on well-understood and familiar programming patterns and reasoning 
techniques. Besides supporting an elegant functional programming style for quan- 
tum computations, the quantum arrow calculus allows reasoning about general 
or mixed quantum computations. This is the first work proposing reasoning 
about mixed quantum computations. The equations of the arrow calculus plus 
the equations of the monadic calculus provide indeed a powerful mechanism to 
make proofs about quantum programs. In [12] we have proposed very similar 
reasoning techniques, however for pure quantum programs. Also, in p^3j the au- 
thor presents a quantum lambda calculus based on linear logic, but just for pure 
quantum computations. 
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A Simply- Typed Lambda Calculus 

The simply-typed lambda calculus with the type of booleans, and with let and 
if is shown in Figure 6. Let A, B, C range over types, L, M, N range over terms, 
and r, A range over environments. A type judgment F \- M : A indicates that 
in environment F term M has type A. As presented in the arrow calculus [3, 
we are using a Curry formulation, eliding types from terms. 



Fig. 6. Simply-typed Lambda Calculus 



Syntax 

Types A, B, C ::= Bool \ A x B \ A B 

Terms L, M, N :■- x \ True | False | (M, iV) | fst i | snd L \ Xx.N \ L M 

let X = M in iV I if L then M else iV 
Environments r,A ::= Xi : Ai, . . . ,Xn : An 
Types 

{x:A)er 



I- False : Bool h True : Bool x ■ A 

r\- M -.A r^N -.B rh L:AxB L:AxB 



r\- {M,N) : Ax B r h fst L : A rh snd L : B 

r,x: Ah N : B F h L : A ^ B F h M : A 



r h Xx.N : A^ B r\- L M : B 

r\- M -.A rx:AhN:B B \- L : Bool Bh M,N : B 



r h let X = Af in TV : B rh if L then M else TV : B 

Laws 

fst {M,N) = M 

(/?!) snd(Tkf,Ar) =TV 

(ry^) (fstL,sndL) =L 

(/?-*) {\x.N)M = N[x ■- M] 

(rj^) \x.{L x) = L 

(let) let a; = TW" in TV = N[x ■- M] 

{/3i) if True then M else TV = M 
(/Ja) if False then A/ else TV = TV 



B Monadic Calculus 



The simply-typed lambda calculus presented in Appendix[A]is the foundation of 
purely functional programming languages. In this section we show the monadic 
calculus [7], which also models monadic effects. A monad is represented using 
a type constructor for computations m and two functions: return :: a ^ m a 
and :>=:: m a ^ {a ^ m b) ^ m b. The operation (pronounced "bind") 
specifies how to sequence computations and return specifies how to lift values 
to computations. From a programming perspective, a monad is a construct to 
structure computations, in a functional environment, in terms of values and 
sequence of computations using those values. 

The monadic calculus extends the simply-typed lambda calculus with the 
constructs in Figure 7. Unit and bind satisfy left unit, right unit, and associa- 
tivity laws, (left), (right), and (assoc). 



Fig. 7. Monadic Calculus 
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MonadPlus Laws 

mzero + a = a 
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a + (6 + c) ={a + h) + c 

let X = mzero in T = mzero 

let X = (M + TV) in T = (let x = M in T) + (let a; = TV in T) 



Beyond the three monad laws discussed above, some monads obey the - 
MonadPlus laws. The MonadPlus interface provides two primitives, mzero and 
+ (called mplus), for expressing choices. The command + introduces a choice 
junction, and mzero denotes failure. 

The precise set of laws that a MonadPlus implementation should satisfy is 
not agreed upon [H], but in |15j is presented a reasonable agreement on the 
laws. We use in Figure 7 the laws introduced by p^ . 

The intuition behind these laws is that MonadPlus is a disjunction of goals 
and is a conjunction of goals. The conjunction evaluates the goals from left- 
to-right and is not symmetric. 

C General Quantum Computations 

Quantum computation, as its classical counterpart, can be seen as processing 
of information using quantum systems. Its basic idea is to encode data using 
quantum bits (qubits). In quantum theory, considering a closed quantum system, 
the qubit is a unit vector living in a complex inner product vector space know as 
Hubert space [1] . We call such a vector a ket (from Dirac 's notation) and denote 
it by \v) ( where v stands for elements of an orthonormal basis), a column vector. 
Differently from the classical bit, the qubit can be in a superposition of the two 
basic states written as q;|0) + or 



a 

/3 

with |ap + |/?p = 1. Intuitively, one can think that a qubit can exist as a 
0, a 1, or simultaneously as both and 1, with numerical coefficient (i.e., the 
probability amplitudes a and (3) which determines the probability of each state. 
The quantum superposition phenomena is responsible for the so called "quantum 
parallelism." 

Operations acting on those isolated or pure quantum states are linear op- 
erations, more specifically unitary matrices S. A matrix A is called unitary if 
S*S = /, where S* is the adjoint of S, and / is the identity. Essentially, those uni- 
tary transformations act on the quantum states by changing their probability 
amplitudes, without loss of information (i.e., they are reversible). The appli- 
cation of a unitary transformation to a state vector is given by usual matrix 
multiplication. 

Unfortunately in this model of quantum computing, it is difficult or impos- 
sible to deal formally with another class of quantum effects, including measure- 
ments, decoherence, or noise. 

Measurements are critical to some quantum algorithms, as they are the only 
way to extract classical information from quantum states. 

A measurement operation projects a quantum state like a\0) + (3\1) onto 
the basis |0),|1). The outcome of the measurement is not deterministic and it 
is given by the probability amplitude, i.e., the probability that the state after 
the measurement is |0) is |ap and the probability that the state is |1) is If 
the value of the qubit is initially unknown, than there is no way to determine a 
and (3 with that single measurement, as the measurement may disturb the state. 
But, after the measurement, the qubit is in a known state; either |0) or 
In fact, the situation is even more complicated: measuring part of a quantum 
state collapses not only the measured part but any other part of the global state 
with which it is entangled. In an entangled state, two or more qubits have to 
be described with reference to each other, even though the individuals may be 
spatially separated El- 
There are several ways to deal with measurements in quantum computing, 
as summarized in our previous work [5]. To deal formally and elegantly with 
measurements, the state of the computation is represented using a density matrix 
and the operations are represented using superoperators [16]. Using these notions, 
the projections necessary to express measurements become expressible within the 
model. 

Intuitively, density matrices can be understood as a statistical perspective of 
the state vector. In the density matrix formalism, a quantum state that used to 
be modeled by a vector \v) is now modeled by its outer product where 
{v\ is the row vector representing the adjoint (or dual) of \v). For instance, the 
state of a quantum bit \v) = "Tf I*-*) "73 1"*^^ represented by the density matrix: 



* For more detailed explanation about entangled, see [l]. 
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Note that the main diagonal shows the classical probability distribution of basic 
quantum states, that is, these state has ^ of probability to be |0) and ^ of 
probability to be |1). 

However, the appeal of density matrices is that they can represent states 
other than the pure ones above. In particular if we perform a measurement on 
the state represented above, we should get |0) with probability 1/2 or |1) with 
probability 1/2. This information, which cannot be expressed using vectors, can 
be represented by the following density matrix: 

1/2 0\ f \ _ f 1/2 

oy + i^o i/2y 1^ 1/2 

Such a density matrix represents a mixed state which corresponds to the 
sum (and then normalization) of the density matrices for the two results of the 
observation. 

The two kinds of quantum operations, namely unitary transformation and 
measurement, can both be expressed with respect to density matrices pTj. Those 
operations now mapping density matrices to density matrices are called super- 
operators. A unitary transformation S maps a pure quantum state \u) to S\u). 
Thus, it maps a pure density matrix \u){u\ to S\u){u\S*. Moreover, a unitary 
transformation extends linearly to mixed states, and thus, it takes any mixed 
density matrix A to SAS*. 

As one can observe in the resulting matrix above, to execute a measurement 
corresponds to setting a certain region of the input density matrix to zero. 



D Definition of Measurement Operations 

In this section we present the denotations of the programs for measurements, trl 
and meas, added to the quantum arrow calculus. 



trL :: Super {A,B) B 

trL((ai, 5i), (02, 62)) = if ai == 02 then return{bi, 62) else mzero 
meas :: Super A {A, A) 

meas(ai, 02) = if == 0-2 then return{{ai, ai), (ai, ai)) else mzero 

We consider projective measurements which are described by a set of projections 
onto mutually orthogonal subspaces. This kind of measurement returns a classi- 
cal value and a post-measurement state of the quantum system. The operation 
meas is defined in such a way that it can encompass both results. Using the 
fact that a classical value m can be represented by the density matrix \m){m\ 
the superoperator meas returns the output of the measurement attached to the 
post-measurement state. 



